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Prufungsantrag gem. § 44 PatG ist gestellt 

(S) Verfahren zur Konturierung von Korperstrukturen in CT-Datensatzen 

(§7) Bei einem Verfahren zur Konturierung von Korperstruk- 
turen in CT-Datensatzen werden zunachst Stutzstellen ge- 
sucht, dann anhand der gesuchten Stutzstellen eine erste 
Modellkontur ermittelt sowie die mit der Modellkontur 
am schlechtesten vertrayliche Stutzstellen entfernt und 
anhand der verbleibenden Stutzstellen eine neue Modell- 
kontur ermittelt, bis die Abweichung der am schlechte- 
sten vertraglichen Stutzstellen einen Schwellwert oder 
die Zahl der verbleibenden Stutzstellen einen Abbruch- 
wert unterschreitet. 
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Beschreibung 

Die Erfindung hetriffi ein Verfahren zur Konturierung von Kdrperstrukturen in CT-Datensatzen, bei welchem eine 
Modellkontur an die tatsachliche Konlur der Korperslruktur angepaBt wird. 

Insbesondere bei der Planting von Tumorbestrahlungen ist die Konturierung von Risikoorganen, wie beispielsweise 
der Augen, von enlscheidender Bedeulung. Diese schaffl die Voraussetzung, da8 die Beslrahlung optimal gewahlt wer- 
den kann, ohne die Risikoorgane niehr als unbedingt nolig zu belasten. 

Bisher muBten die Risikoorgane manuell mittels graphischer Werkzeuge konturiert werden. Allein fur die Augen sind 
hierrur pro Patient in etwa 6 min. zu veranschlagen. 

Es ist Aufgabe vorliegender Erfindung eine automaiische oder semiautornatische Konturierung von Korperstrukluren 
in CT-Datensatzen mit ausreichender bzw. beherrschbarer Zuverlassigkeit zu ermoglichen. 

Als Losung schlagt die Erfindung ein Verfahren zur Konturierung von Korperstrukluren in CT-Datensatzen vor bei 
welchem eine Modellkontur an die tatsachliche Kontur der Korperstrukluren angepaBl wird, wobei zunachst Stut7stellen 
gesucht werden, anhand der gesuchten Stutzstellen eine erste Modellkontur ermittelt wird sowie die mit der Modellkon- 
tur am schlechtesren vertragliche Stutzstelle entfemt und anhand der verbleibenden Stutzstellen eine neue Modellkontur 
ermittelt wird, bis die Abweichung der am schlechtesten vertraglichen Stutzstelle einen Schwellwert oder die Zahl der 
verbleibenden Stutzstellen einen Abbruchwert unterschreitet. 

Das auf diese Weise bereitgestellte Verfahren erweist sich als auBerst robust gegenuber fehlerhaften Stutzstellen da 
nicht alle ermitlelten Stutzstellen sondern nur Stutzstellen innerhalb einer bestimmten Fehlerbandbreite zur Besrimmunc 
20 der Modellkontur genutzt werden. 

In der Praxis hat sich gezeigt. daB bei spiels weise bei einem Auge bei geeigneter Wahl des Schwellwerles in den aller- 
meisten Fallen mehr als 30 Stutzstellen die Eliminationsroutine passierten. Selbst bei einem Anteil fehlerhafter Stutzstel- 

lC Y,°™ 35% k5nnen dlCSe erkannt und entfemt werden. Der ubliche Anteil fehlerhafter Stutzstellen liegt zwischen 6% 
und 12%. 

25 Es versteht sich, daB statt eines CT-Datensatzes auch andere Bildinformationen enthaltende Datensatze Verwendune 
fmden konnen, bei denen Slrukluren durch eine Modellkontur dargestelll werden sollen. 

Das Verfahren cnnoglichl cs ? cine ansich als Satz von Bildpunktcn vorlicgcndc Datcnincngc, wclchc die Kontur cincr 
Korperslruktur darstellt, auf eine nur wenige Parameter umfassende Modellkontur zu reduzieren. 

Beispielsweise kann eine bestimmte Zahl von Stutzstellen einer Augenkontur auf die Mittelpunktslage sowie die 
Halbachsen einer Ellipse reduziert werden. Wie unmitlelbar ersichtlich, bedingt die Zahl der fur die jeweilige Modell- 
kontur benotigien Parameter (z. B. Mittelpunktslage und Halbachse) den Abbruchwert, ab welchen'das Verfahren als 
miBlungen abgebrochen wird. Sicherheitshalber kann dieser Abbruchwert aber auch hoher gewahlt werden 

Die Genauigkeit des Verfahren laBt sich dadurch erhohen, daB der CT-Datensatz vor der Stutzstellensuche durch einen 
Filter vorbearbeitet wird. 

Insbesondere ist es moglich, als Filter eine von der Hesse-Matrix abhangige Funktion auf den CT-Datensatz anzu wen- 
den^ Ein derartiger Filter kann derartig gewahlt werden, daB er den CT-Datensatz entlang der Eigenrichtung des betrags- 
maBig groBeren Eigenwertes der Hesse-Matrix scharft und senkrecht dazu entlang der Eigenrichtung des betragsmaBig 
kleineren Eigenwertes glattet. Somit hebt der Filter linienarlige Strukturen hervor. Beispielsweise kann der Filter die bei- 
den unabhangigen SO(2)-Invarianten der Hesse-Matrix Tr H und Tr H 2 umfassen. Diese SO(2Hnvarianten sind gegen- 
uber Raumdrehungen invanant, so daB diese unabhangig von einer konkreten Orientierung angewandt werden konnen 
Auf diese Weise wird sichergestellt, daB die durch den Filter hervorgehobenen Strukturen unabhangig davon sind wie 
die Lage des Korpers im CT-Datensatz ist. Bei diesen Funktionen ist eine Berechnung der Eigenvektoren und Eigenwerte 
der Hesse-Matnx nicht notwendig, so daB der Filter verhaltnismaBig schnell auf den CT-Datensatz angewandt werden 
kann. 

Die Stutzstellensuche erleichtert sich, wenn sie von einen Ort im Inneren der zu konturierenden Korperslruktur aus- 
geht. Bei einem derartigen Ausgangsort ist in der Regel der Rand der Korperslruktur als von innen gesehen erste Struktur 
deiekuerbar Einzelne dazwischen liegende Strukturen oder Schattierungen konnen durch geeignete Auswerteverfahren 
wie beispielsweise das erfindungsgemaBe Konturierungs verfahren, eliminiert werden. Andererseits ist eine derartige 
Wahl eines Ortes, der als Ausgangspunkt genutzt werden kann, auch unabhangig von den Merkmalen des Konturierungs- 
50 vertahrens t Lir die Bildauswertung vorteilhaft . 

Die Oriswahl kann, insbesondere bei semiautomatischen Verfahren, durch den Nutzer erfolgen, indem dieser einen 
enisprechenden Ort im Inneren der zu konturierenden Korperslruktur markiert. Es ist andererseits auch denkbar daB ein 
derartiger Ort durch geeignete zusatzliche Verfahrensschritte ermittelt wird. 

Hierzu kann zum Beispiel in der Umgebung eines jeden Punktes, beispielsweise in einer Umgebung von 7 x 7 Bild- 
punkten, des CT-Datensatzes an die gefilterten Hounsfleld-Werte eine Funktion F angepasst werden, die die Gestalt eines 
gekrummten Kammes hat, wobei die Krummung der Kammlinie konslant sein soli. Als Kammlinien-Funktion G kann 
hierbei eine behebige Funktion, die achsensymmetrisch ist und von Null ausgehend abfallt, dienen. Die Kammlinien- 
Tunktion G kann beispielsweise eine Glockenkurve, wie z. B. eine GauB-Funktion sein. Insofern ist die anzupassende 
Funktion F in ihrem Krummungsradius sowie dem Mittelpunkt der Krummung anpassbar. Daruber hinaus konnen noch 
konslante Fit-Parameter hinzugefugt werden. Auf diese Weise kann fur jede Umgebung eines jeden Punktes eine Kamm- 
hme enriittelt werden, die ein Kreissegment mit einem Mittelpunkt darsteUt. Indem die Mittelpunktskoordinaten fur je- 
den Punkt des CT-Datensatzes berechnet werden, laBt sich eine Haufigkeitsverteilung der Miltelpunkte ermitteln Urn 
eine Verwechslung mil anderen Strukturen des CT-Datensatzes zu vermeiden, werden in die Verteilung lediglich dieie- 
nigen MitteJpunkte aurgenommen, bei denen der jeweilige Krummungsradius in etwa dem Radius der Korperstrukluren 
entspnent, beispielsweise zwischen einem minimalen und einem maximalen Augenradius liegt Die auf diese Weise ge- 
wonnene Verteilung der Mittelpunktskoordinaten enthalt zwei dominante Peaks, die beispielsweise den Mittelpunkten 
der beiden Augen entsprechen. Insofern braucht die Verteilung lediglich geglattet, der Untergrund entfemt und die Lage 
der Peaks ermittelt werden, urn einen Ort in der Korperslruktur zu ermitteln. 
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Zur Stulzstellensuche konnen von einem im Innercn der Korpersiruklur ausgewahlten On ausgehend entlang eindi- 
mensionaler Wege lokale Maxima des CT-Datensatzes enuiltelt werden, bei denen die Differenz zwischen einem voraus- 
gegangenen Minimum und deni jeweiligen Maximum zwischen einem unleren Siuizslellenschwellwert und einem obe- 
ren Siiiizstellenschwellwert liegt. Hierbei bezeichnei der BegrifT eindimensionaler Weg einen von dem gewiihlien Ort 
ausgehenden Weg durch das zweidimensionale Bild des CT-Datensatz.es beliebiger Art, der keine Uherkreuznngen mil 5 
sich selbst aufweist. 

Das Verfahren gestaliet sich vorteilhaft, wenn Strecken verwendet werden, so daB schleifende Schnille mil der Kontur 
vermieden werden. 

Die Differenz zwischen Minimum und Maximum erweist sich als ein verhallnismaBig einfach umzusetzendes und 
dennoch zuverlassiges Kriteriuni fur die Stutzsteliensuche. Hierbei konnen der untere Stutzsiellenschwellwert und der to 
obereStutzstellenschwellwert verhaltnismaBig unkritisch als Verfahren sparameier gewahlt werden, so daB eine derartige 
Stutzsteliensuche auch unabhangig von den iibrigen Merkmalen des Konturierungsverfahrens vorteilhaft ist. Insbeson- 
dere im Zusammenspiel mit dem erfindungsgemaBen Konturierungs verfahren erweist sich diese Stutzsteliensuche als in 
der Wahl der Verfahrensparameler unkritisch, da fehlerhafte Stutzstellen verhaltnismaBig zuverlassig eliminiert werden 
konnen. 15 

Anstelle einer Differenz zwischen Minimum und Maximum kann auch eine Steigung zwischen Minimum und Maxi- 
mum als Kriteriuni fur die Stutzsteliensuche dienen. Hierbei solite die Steigung ebenfalis zwischen zwei Stutzstellen- 
schwellwerten liegen und als Stutzstellenkrirerium genurzt werden. Auch hier ergeben sich die vorbeschriebenen Vor- 
teile. Ebenso kann ein absoluter Schwellwert. angegeben werden, der iiberschritten werden muB, wobei dann das nachste 
Maximum als maBgebliche GroBe verarbeitet wird. 20 

Um zu vermeiden, daB die Konturierung der Korperstrukiuren zu klein erfolgt und hierdurch unbeabsichtigt Schaden 
bedingt werden, kann nach der vorbeschriebenen groben Stutzsteliensuche als Feinjusl.ierung bei spiels Wei se ein lokales 
Maximum und ein nachstliegendes, von der Korperstruktur ausgehend entlang eines eindirnensionalen Weges weiter au- 
Ben liegendes lokales Minimum des CT-Datensatzes gewahlt werden. Bei einer derartigen Anordnung kann die Stutz- 
stelle durch eine geeignete bzw. durch eine uber einen Parameter bestimmbare Interpolation zwischen dem Maximum 25 
und dem nachst auBenliegenden Minimum definiert werden. 

Es vcrstcht sich, daB bei den vorbeschriebenen Anwcndungcn auf den CT-Datcnsatz nicht nur ein CT-Datcnsatz von 
Hounsfield-Werten Verwendung finden kann, sondern daB hier auch andere Bilddaten bzw. auch durch eine Filter vorbe- 
handelte CT-Daten oder Bilddaten verarbeitet werden konnen. 

Eine zuverlassige Stutzsteliensuche laBt sich beispielsweise dadurch gewahrleislen, daB hierfur eine Kurve verwendet 30 
wird, die in einem Bereich zwischen dem iokalen Maximum und dem nach auBen hin nachstliegenden lokalen Minimum 
zwischen einem Polynom 3. Grades, welches in einem gewahlten Inlervall den CT-Datensatz approximiert und dessen 
quadratischer Naherung in diesem Intervall liegt. 

Hierbei kann das Intervall beispielsweise die Halftedes Bereiches zwischen den aufeinanderfolgenden, das lokale Ma- 
ximum einrahmenden Minima ausmachen. 35 

Vorzugsweise ist die Kurve eine Interpolation zwischen dem Polynom 3. Grades und dessen quadratischer Naherung 
um das Maximum des Polynoms 3. Grades. Diese Kurve kombiniert die Vorteile einer Stutzstellenbestimmung mittels 
Naherungen 2. und 3. Grades an den CT-Datensatz, besitzt jedoch die Nachteile dieser Approximationen nicht. So ist die 
Lage der Stutzstelle wie bei der Verwendung einer Polynomapproximation 2. Grades immer definiert, was bei einer Po- 
lynomapproximation 3. Grades nicht zwingend der Fall sein muB. Daruber hinaus berucksichtigt diese Kurve eine 40 
Schiefe des CT-Datensatzes um das lokale Maximum ahnlich wie eine Polynomapproximation 3. Grades. Dieses konnte 
eine Approximation 2. Grades nicht leisten. 

Die Verwendung einer derartigen Kurve ist auch unabhangig von den iibrigen Merkmalen des Konturierungsverfah- 
rens fiir eine Stutzsteliensuche vorteilhaft anwendbar. Insbesondere kann eine derartige Stutzsteliensuche auch mit ande- 
ren als dem erfindungsgemaBen Verfahren fiir eine Modellkonturermittlung kombiniert werden. 45 

Vor der Ennittlung der ersten Modellkontur konnen Teilstrukturen der Korperstruktur errnittelt und entsprechende 
Stutzstellen markiert bzw. selektiert werden. Auf diese Weise beeinflussen diese Stutzstellen die nachfolgende Modell- 
konturenniulung nicht nachteilig. Beispielsweise konnen bei einem Auge vorab die Cornea (Homhaut vor der Linse) 
bzw. der Sehnerv segmentiert werden. Fiir die Ennittlung der Teilstruktur konnen beispielsweise die vorbeschriebenen 
Verfahren Anwendung finden, wie sie anhand der Konturierung der Gesamt korperstruktur beschrieben sind. 50 

Es versteht sich, daB das erfindungsgemaBe Verfahren im Prinzip in beliebigen Dimensionen, insbesondere 2- und 3- 
diniensional, Anwendung finden kann. Hierbei ist die 3-dimensionale Anwendung auch unabhangig von den iibrigen 
Merkmalen des Konturierungsverfahrens vorteilhaft, da hierbei von einem Ort im Inneren der Korperstrukiuren ausge- 
hend bei jedem Schnitt ein Schleifen entlang der Konlur vermieden werden kann. 

Weitere Vorteile, Ziele und Eigenschaften vorliegender Erfindung werden anhand der Beschreibung anliegender 55 
Zeichnung dargestellt, in welcher beispielhaft ein errindungsgemaBes Konturiervcrfahren dargestellt ist. 

In der Zeichnung zeigen: 

Fig. 1 eine schematische Darstellung des Verfahrensflusses bei einem erfindungsgemaBen Konturierungs verfahren, 
Fig. 2 eine schematische Darstellung einer Lageparameterintervallbestinimung bei dem Verfahren nach Fig. 1, 
Fig. 3 eine Prinzipskizze zur Grobhestimmung der Stutzstellen lage bei der Stutzsteliensuche, 60 
Fig. 4 eine vergroBerte Darstellung der Skizze nach Fig. 3 zur Darstellung einer Feinjustierung bei der Stutzsteliensu- 
che, 

Fig. 5 eine beispielhafte Darstellung eines CT-Datensatzes uber einem eindirnensionalen Weg, 

Fig. 6 gefundene Stutzstellen und emiittelte Modellkontur bei einem ersten Paar aus unterem und oberen Stutzstellen- 
schwellwert und 65 

Fig. 7 gefundene Stutzstellen und ermittelte Modellkontur bei einem zweiten Paar aus unterem und oberen Sliitzstel- 
lenschwellwert. 

Das in den Figuren dargestellte Konturierungs verfahren dient der Anpassung einer Ellipsenform an eine tatsachliche 
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Kontur eines Auges. Hierzu werden, wie in Fig, 1 dargestellt, die CT-Daten zunachsTriner Vbrbearbeitung durch einen 
Filler untcrzogen. In einem zweilen Arbeit sschritt werden zunachst im Rahmen einer Grobsuche Stutzstellen auf dem je- 
wcihgcn Augcnrand gesuchl. Im AnschluB hieran wird anhand der gefundenen Slutzstellen die Linse gesuchl und die 
Cornea neu segmentiert. Bevor anhand der gesuchten Stutzstellen eine erste Modellkontur ermittelt wird, uni fehlerhafte 
Stulzslellen zu climimeren, findel eine Feinjustierung der Stulzslellen statt. Im AnschluB hieran wird wie bereits vorste- 
hend angedculcl, eine Modellkontur ermittelt, wobei fehlerhafte Slutzstellen eiiminierl werden Nach einer abschlieBen- 
den Llhpsenanpassung wird eine auf diese Weise als Modellkontur gewonnenen Ellipse als Kontur ausgegeben 

Bet dem ycrwcndclen Filler wird eine Funktion aus den beiden unabhangigen SO(2)-Invarianten der Hesse-Matrix Tr 
H und 1 r 1 \~ 

h'(x,y) = h(x,y ) a /2Tr(TT(x^p) - Tr-II(x,y)sign(TrII(x,y)); 0<a< 1 

auf die Hotinsheld-Werlc der zweidimensionalen CT-Daten angewandt. 

Dieses isi beispielhafi fur einen Schnitt entlang eines eindimensionalen Weges durch einen CT-Datensatz in Fig 5 dar- 
gestellt. Hierbei bezcichnei h die Hounsfield-Werte, wahrend h' die Werte nach der Filterung darstellt Der Filter hebt in 
dem zweidimensionalen Daiensatz iinienartige Strukturen hervor, laBt jedoch einzeln stehende Bildpunkte, die sich von 
der Umgebung abheben. unvcrandert. Auf diese Weise werden in dem CT-Datensatz vorhandene, Iinienartige Konturen 
verscharti. Wie unmmelbar aus Fig. 5 ersichtlich, werden somit Anstiege uberhoht. In einem derartigen Schnitt ist ein 
Olatten entlang der limcnarligcn Strukturen nicht darstellbar. Der Filter kann itteriert werden. 

In dem nachsten Verlahrcnsschritt werden Stutzstellen grob ermittelt. Hierbei wird zunachst ein Ort im Inneren des 
Auges gewahll. Dieses geschieht manuell, so daB es sich bei den hier beschriebenen Verfahren urn ein semiautomatisches 
Kontunerverfahren handell. hs hai sich gezeigt, daB das vorliegende Konturierverfahren verhaltnismaBig unempfindlich 
gegenuber einer Ortswahl in dieser Hinsicht ist, so daB der Ort lediglich innerhalb einer Ellipse mit halb so groBen Halb- 
achsen wie die des Auges gewahll werden sollte, urn ausreichend gute Ergebnisse zu erzielen. 

In einer weiteren Ausgest all ung dieses Verfahrens konnten der Ort im Inneren des Auges bzw. die Orte im Inneren der 
Augen automatiscli cnimiell werden. Hierzu wird an die Umgebung, beispielsweise die Umgebung von 7 x 7 Bildpunk- 
tcn cines jeden Punkies ( x l . y l; ) des CT-Datcnsatzcs cine Funktion F an die gcfiltcrtcn Hounsficld- Werte h' angepasst 
die die Gestalt eines gekm mmien Kammes hat, wobei die Krummung der Kammlinie konstant ist. Die Funktion F hat die 
Form: 



F(x,y) = a G 



l 



+ c 



wobei a, c, x z y 2 und r Parameter sind, die angepasst werden. Hierbei sind (x z , yz ) die Koordinaten des Krummungsmit- 
telpunktes und r der Krummungsradius. Die Kammlinienfunktion G veist hierbei folgende Eigenschaften auf- G (t) = G 
(-0, G*(0 > 0, fur t < 0. Bei vorliegendem Ausfuhrungsbeispiel wird eine auf den Maximalwerr 1 festgelegte GauB'sche 
Glockenkurve gewahll. Fur jedc Umgebung eines jeden Punktes stellt somit die Kammlinie ein Kreissegment dar, dessen 
Mittelpunkt jeweils bei (x/, y y ) liegt. 

Indem die Mittelpunklskoordinaien (x z , y z ) fur samtliche Punkle (x D , yu ) des CT-Datensatzies berechnet werden kann 
erne Haufigkeitsverteilung der Mittelpunkte ermittelt werden. Urn lediglich Kammlinien zu erhalten, die zu den Aueen 
gehoren, gehen lediglich Krummungsmittelpunkte (x z ( Xu , y v \ y z ( Xu , y 0 )) in die Verteilung ein, deren Kilimmunesra- 
dius r zwischen einem gewahllen Minimum und einem gewahllen Maximum liegt. 

Die so ermittelte Verteilung der Mittelpunktskoordinaten weist zwei dominante Peaks auf, die den Mittelpunkten der 
beiden Augen entsprechen. Die Verteilung wird anschlieBend geglattet, der Untergrund entfernt und die Lage der Peaks 
ermittelt. Diese Peaks werden als Startwerte bzw. als Ausgangsort fur die Stutzschwellensuche gewahlt 

Von dem Ausgangsort ausgehend werden eindimensionale Schnitte radial nach auBen durch den CT-Datensatz eeleet 
Derartige Schnitte sind beispielhafi in Fig. 5 und in den Fig. 3 und 4 skizziert. Als nachstes werden von diesem Ort der 
mil 0 bezeichnet ist, ausgehend lokale Maxima x L ermittelt, bei denen die Differenz Ah zwischen dem vorausgegangenen 
Minimum x 0 und dem Maximum x { zwischen einem unteren Stutzslellenschwellwert und einem oberen Stutzstellen- 
schwellwert liegt. Dieses isi in Fig. 3 dargestellt. 

Nachdem auf diese Weise Si iitzstellen gefunden wurden, werden die Teilstmkturen, namlich die Cornea des Auges er- 
mittelt Dieses ertolgt dadurch, daB von dem Ausgangsort in dem Auge ausgehend eindimensionale Schnitte radial nach 
auBen durch den CT-Datensatz gelegt werden, wobei nicht nur das diesem Ausgangsort zunachst liegende lokale Maxi- 
mum sondem auch ein hierauf folgendes lokales Maximum, welches derselben Bedingung geniigt, gewahll wird Diese 
Schnitte werden jedoch nur m einem nach vome weisenden Winkelsegment von 120° gewahlt, indem eine Linse zu er- 
warten ist. Dieses ist moglich, da die Lage eines Kopfes in dem CT-Datensatz bekannt ist. Andererseits ist es denkbar 
durch weitere Vertahrensschritle, ahnlich der bei der Ermittlung der Augenmittelpunkte verwendeten, die Laee dieses 
Winkelsegmentes zu best im men. 

Im AnschluB hieran wird von der im Uhrzeigersinn liegenden Seite des Winkelsegments ausgehend die let/te Stiitz- 
stelle So gewahlt, ab welcher entgegen dem Uhrzeigersinn gesehen zwei Altemativen S, und Sr bestehen Diese beiden 
Alternativen werden jeweils uber eine Linie mil S 0 verbunden und jeweils die kleinsten Hounsfield-Werte entlang dieser 
Linien ermittelt. Als Bestandteil der zu selektierenden Kontur wird diejenige Stutzstelle gewahlt, bei welcher der kleinste 
Hounsheld-Wert h' groBer ist. Von der auf diese Weise ausgewahlten StutzsteUe ausgehend, wird itterativ das entgegen 
dem Uhrzeigersinn nachstliegende Stutzstellenpaar gewahlt und von diesen in derselben Weise eine Stutzstelle als Be- 
standteil der Kontur gewahlt. 

Dieser Vorgang wird solange wiederholt, bis das gesamte Winkelsegment von 120° durchlaufen ist bzw. bis lediglich 
nur noch eine Alternative vorliegt. & 
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Aufdiesc Wcise wird eine linienarlige Konlur emii licit, die moglichsl kontinuierlich von dem bereits bcslinimten Au- 
gcnrand, welcher sich durch lediglich eine vorliegende Alternative kennzeichnet, ausgehend weiterverlaufl. 

1st die aufdiesc Weise ermitteiie Alternative die auBenliegende, so handelt es sich bei dieser Alternative um die Cor- 
nea. Die innenliegendc Alternative ist dann die Lin se. 

1st riicsc Alternative die inncnliegende, so gehort die auBere Alternative zum Lid odcr zu einer anderen, nicht zum 5 
Auge gehorenden Struktur. 

Aufdiesc Weise ergibl sich eine Menge an Stutzstellenwerten, die den Augenrand sowie die Cornea umfassen und die 
in cincin nachsien Verfahrcnsschritt. feinjusliert werden. Hierzu wird um das lokale Maximum x\ in einem Intervall, wel- 
ches die I lalt'icdcs Bereiches zwischen dent vorausgegangenen Minimum Xq und dent nachlolgenden Minimum x 2 aus- 
muchi an IV ein Polynoin 3. Grades io 

P3(x)=ax* + bx 2 + cx + d 

odcr in i Koordinatensystem des lokalen Maximums x M ^ Xj 
PM\ x M ) = a«(x - x M ) 3 + a 2 (x - x M ) 2 + cto 
mil der quadrat ischen Naherung 
P2(x \ M ) = rx>(x x M ) 2 +.ao 



angepaBl (siche Fig. 4). Zur Berechnung der entsprechenden Stutzstellen werden nunmehr die Schnittpunkte dieser Po- 
lynomapproximalion P3 mil einer konstanten Funktion C = h'(xo) + IlAh in Abhangigkeit vom Koeffizient <x 3 des kubi- 
schen Terms dcr Polynomapproximation P3 berechnet. Von diesen Schnittpunkten wird der auf der AuBenseite des Au- 
ges am wcilcstcn inncn liegende ausgewahlt. Dessen Lage wird anschlieBend im Koeffizienten a 3 Taylor-en twickelt. Zur 
Bcsiimmurig dcr Siut/.sicllenlage l St wird das Taylorpolynom 1. Grades in 0,3 
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verwendct. ITicrbci lie fern die Tupel [1 = lstOl), h' = h'(x 0 ) + IlAh] mil variables Parameter T| eine Kurve K, welche zwi- 35 
schen dcr Polynomapproximation 3. Grades P3 und deren quadratischer Naherung P2 verlauft. 

In Abhangigkeit von dem Parameter T| kann somil eine Stutzstelle bestimmt werden. Hierbei kann uber den Parameter 
T[ cingcsicllt werden, wic cng die AuBenkonlur am Auge sitzen soli. Bei TJ = 1 sitzt die AuBenkontur auf dem Ring der lo- 
kalen Maxima, wahrend sic bei T| = 0 wesentlich weiter auBen in der Nahe des nachstauBeren Minimums liegt (sieheFig. 
4). In der Praxis kann hierdurch eine Stutzstelle um 2 bis 3 mm angepasst werden, wobei der behandelnde Arzt uber TJ 40 
den CJrad dcr Anpassung vorgeben kann. Ublicherweise wird von Arzten ein Wert von T| =s 0,5 gewahlt. 

In einem wcitcrcn Vcrfahrensschritt. wird dann an die verbleibenden bzw. ermittelten Stutzstellen eine Ellipse ange- 
paBt. Im AnschluB hicran wird die Abweichung der schlechteslen Stutzstelle von der Ellipse bestimmt und so lange die 
jcwcils schlechtesle Stiitzstelle entfernt und daraufhin eine neue Ellipse angepaBt, bis der Fehler der schlechtesten Stiitz- 
stelle cincn bcstimintcn, durch die Auflosungsgrenze der CT-Daten festgelegten Wert unterschreitet. 45 

Das Verfahren bricht spateslens ab, wenn nur noch 4 Stutzstellen ubrig sind. Die Konturfindung gilt jedoch bereits 
dann als miBlungcn, wenn von einer anfanglichen Zahl Stutzstellen weniger als eine kritische Zahl verbleiben. So kann 
beispiclswcisc bei anfanglichen 36 Stutzstellen eine kritische Zahl von 8 gewahlt werden. Bei Testlaufen ist dieses je- 
doch nie aufgeirctcn. In den allermeisten Fallen passicren niehr als 88% der Stutzstellen die Eliminationsroutine, Die 
Routine fiir fchlcrhaflcn St iitzslellen ist robust und elTcktiv. Sie kann selbst bei einem Anteil fehlerhafter Stutzstellen von 50 
35% dicsc crkennen und entfernen. Der ubliche Anteil fehlerhafter Stutzstellen liegt zwischen 6% und 12%, Es versteht 
sich, daB .'inch cine andere Zahl der St iitzslellen gewahlt werden kann. Das Verfahren ist von der genauen Zahl der Stiitz- 
stellen in wcitcn Grenzen unabhangig. 

Das Ergchnis dcr Hliminationsrouiine ist in Fig. 6b und 7b beispielhaft dargestellt, wahrend Fig. 6a und 7a als 
schwar/c Kreuzc die verwendcten und als wciBc Kreuzc die eliminierten Stutzstellen zeigen. Hierbei zeigen die Fig. 6a 55 
und 6b eine Vcrfahrcnsfuhrung, bei welcher die Slulzsiellenschwellwerte zur Ermiltlung der lokalen Maxima x t und Mi- 
nima x 0 bzw. x : vcrhaltnismaBig gut eingestellt sind. Fig. 7a und 7b hingegen zeigen ein Ausfuhrungsbeispiel, bei wel- 
chem die Stutzstellcnschwellwerte nicht so gut eingestellt sind. Hierdurch bedingt, weichen die Stutzstellen stark von- 
einander ab. das hci&i es sind vcrhaltnismaBig viele fehlerhafte Stutzslcllen ermitteli worden. Wie jedoch a us Fig. 7b er- 
sichllich, licfert selbst in einem derartigen Fallc das erfindungsgema'Bc Verfahren eine hervorragende Konturierung des C*) 
Augcs. 

Es versteht sich. daB die in dem Ausfuhrungsbeispiel dargcslelllc Verfahrensfolge kleineren Abweichungen untcrzo- 
gen werden kann. So kann beispiclswcisc aus dcr erstnialigcn Segmcntierung des Augenrandes zunachst ein erstes, noch 
relativ grobes Modcll enuiUell werden, bevor eine Segmcntierung von Cornea und Linse vorgenommen wird. Ebenso ist 
es nioglich, bereits nach der ersien groben Siutzstellensuche eine Feinjustierung vorzunehmen, bevor eine Segmentie- 65 
rung von Cornea und Linse vorgenommen wird. 
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Patent anspruche 

1 . Verfahren zur Konturierung von Korperstrukturen in CT-Datensatzen, bei welchem eine Modellkontur an die tat- 
sachlichc Kontur der Korperstruktur angepasst wird, dadurch gekennzeichnet, 

5 - daB zunachst Stutzstellen gesucht werden, 

- daB anhand der gesuchten Stutzstellen eine erste Modellkontur ermittelt wird, 

- daB die mil der Modellkontur am schlechtesten vertragliche Stutzstelle entfernl und anhand der verbleiben- 
den Stutzstellen eine neue Modellkontur ermittelt wird, bis die Abweichung der am schlechtesten vertragli- 
chen Stutzstelle einen Schwellwert oder die Zahl der verbleibenden Stiitzstellen einen Abbruchwert unter- 

10 schreitet. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB der CT-Datensatz vor der Stutzstellensuche durch ei- 
nen Filter vorbearbeitet wird. 

3. Konturierungs verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB zur Stutzstellensuche zunachst 
ein Ort im Inneren der Korperstruktur gewahlt wird. 

15 4. Konturierungsverfahren nach Anspruch 3, dadurch gekennzeichnet, daB zur Stutzstellensuche von diesem Ort 

ausgehend entlang eindimensionaler Wege lokale Maxima des CT-Datensatzes ermittelt werden, bei denen die Dif- 
ferenz zwischen einem vorausgegangenen Minimum und dem jeweiligen Maximum zwischen einem unteren Stutz- 
stellenschwellwert und einem oberen Stutzstellen schwellwert liegr. 

5. Konturierungsverfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, daB zur Stutzstellensuche 
20 ein lokales Maximum und ein nachstliegendes, von dem Inneren der Korperstruktur ausgehend entlang eines eindi- 

mensionalen Weges weiter innen liegendes Minimum des CT-Datentragers gewahlt werden. 

6. Konturierungsverfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, daB zur Stutzstellensuche 
ein lokales Maximum und ein nachstliegendes, von der Korperstruktur ausgehend entlang eines eindimensionalen 
Weges weiter auBen liegendes lokales Minimum des CT-Datensatzes gewahlt werden. 

25 7. Konturierungsverfahren nach Anspruch 6, dadurch gekennzeichnet, daB zur Stutzstellensuche eine Kurve ver- 

wendel wird, die in einem Bereich zwischen dem lokalen Maximum und dem nach auBen hin nachslliegenden lo- 
kalcn Minimum zwischen einem Polynom 3. Grades, welches in cincm gcwahlten Intcrvall den CT-Datcnsarz ap- 
proximiert, und dessen quadratischer Naherung in diesem Intervall liegt. 

8. Konturierungsverfahren nach Anspruch 7, dadurch gekennzeichnet, daB die Kurve eine Interpolation zwischen 
30 dem Polynom 3. Grades und dessen quadratischer Naherung ist. 
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